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Abstract 

The hemisphere soft function is calculated to order a 2 s . This is the first multi-scale soft 
function calculated to two loops. The renormalization scale dependence of the result 
agrees exactly with the prediction from effective field theory. This fixes the unknown 
coefficients of the singular parts of the two- loop thrust and heavy-jet mass distributions. 
There are four such coefficients, for 2 event shapes and 2 color structures, which are 
shown to be in excellent agreement with previous numerical extraction. The asymptotic 
behavior of the soft function has double logs in the CfCa color structure, which agree 
with non- global log calculations, but also has sub- leading single logs for both the CfCa 
and CpTprij color structures. The general form of the soft function is complicated, does 
not factorize in a simple way, and disagrees with the Hoang-Kluth ansatz. The exact 
hemisphere soft function will remove one source of uncertainty on the a s fits from e + e~ 
event shapes. 



1 Introduction 



There has been significant activity in the last few years in the effective field community to 
perform accurate calculations of event shapes for e + e~ colliders. At high energy, the hadronic 
final states in e + e~ collisions are dominated by the formation of jets of particles and are 
described by perturbative QCD. Comparison of theoretic calculations of event shapes with 
the experimentally measured values has lead to some of the most precise measurements of 
the strong coupling constant a s . The NNLO fixed order calculations in [1, 2, 3, 4] allow the 
prediction of many events shapes to order a^. Advances in Soft-Collinear Effective Theory 
(SCET) [5, 6, 7, 8] have allowed for resummation of large logarithmic corrections to thrust 
[9, 10, 11] and heavy jet mass [12] to N 3 LL accuracy and non-perturbative considerations 
were included for thrust in [13]. These results have been used to extract a value of a s that is 
competitive with the world average [14]. 

Dijet event shapes such as thrust and heavy jet mass demonstrate singular behavior when 
calculated perturbatively at fixed order due to the appearance of large logarithmic corrections. 
These large logarithms invalidate a naive expansion in a s and thus need to be resummed to 
provide accurate predictions in the dijet limit. In the dijet limit, there is a clear separation 
between scales. Effective theory techniques rely on a separation between kinematic scales and, 
through renormalization group (RG) evolution, logarithms of the ratio of these scales can be 
resummed. Each of the relevant scales is described by different physics, each of which can 
be calculated using a different theory. The contribution from each can be shown to factorize 
into a hard contribution, due to physics at the center of mass energy Q, a jet function, due to 
physics at the jet scale, and a soft function which describes soft gluon emission. The hard and 
jet functions are known to 2-loops. However, the soft function relevant for thrust or heavy 
jet mass is only partially known beyond 1-loop [15, 12]. In this paper, the perturbative soft 
function is computed analytically to order a 2 .. 

Soft functions have been studied for many years, not just in SCET. These soft functions 
are defined as matrix elements of Wilson lines. For resummation up to the next-to-leading 
logarithmic order (NLL), all that is needed about the soft function is its anomalous dimension. 
This can be extracted either from renormalization-group invariance or from the virtual graphs. 
For example, such calculations have been done for thrust [11], direct photon [16], and dijet 
production [17, 18, 20, 21]. To go beyond NLL, one needs the finite parts of these soft 
functions, which are more difficult to calculate because the real emission graphs are needed, 
and these involve often complicated phase-space cuts. In all cases calculated at 2-loops so 
far, such as Drell-Yan [22, 23] or b — > 57 [24], the real emission graphs only involve one scale. 
Multi-scale soft functions, where different constraints are placed on gluons or quarks going in 
different directions, such as the hemisphere soft function, are likely to play an important role 
in hadron collisions [19, 25]. At order a s , the multiple scales are irrelevant, since only one 
gluon can be emitted. At order a 2 or beyond, there can be real emission graphs depending 
on multiple scales at the same time. It has been suggested [15] that the soft function should 
depend only on logarithms of these scales, such as ln 2 (/c^//c/j). Whether more complicated 
scale-independent terms, such as Li 2 (— ^//cr) + Li 2 (— kji/ki) might appear has been an open 
question. Understanding the form of these soft functions in more detail will be important for 
LHC precision jet physics at NNLL and beyond [25]. 
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The hemisphere soft function S(ki, kn, ll) is the probability to have soft radiation with 
small component k^ going into the left hemisphere and soft radiation with small component 
kfi going into the right hemisphere. More precisely, in e + e~ — > hadron events at center-of-mass 
energy Q, in the limit that all radiation is much softer than Q, the cross section is given by 
matrix elements of Wilson lines. These Wilson lines point in the direction of two back-to-back 
light-like quarks which come from the Born process e + e~ — > qq. Each quark direction defines 
a hemisphere, which we call left and right and denote with the light-like 4-vectors n 11 and 
n M . If the total radiation in the left (right) hemisphere is P£ (-Pr), then S(ki,, k^, /j,) is the 
matrix element squared to have k^ = n - P^ and &r = n- Pr, with all other degrees of freedom 
integrated over. 

The hemisphere soft function is known to have many interesting properties and is con- 
jectured to have others. The factorization theorem for the full hemisphere mass distribution 
implies that the Laplace transform of the soft function should factorize into the form 

s(L u L 2 , Li) = S M (Li)s A1 (L 2 )s / (Li - L 2 ) (1) 

where L\ = Inxip an d L 2 = hxxuLi, with xr and xl the Laplace conjugate variables to ki 
and /cr. The anomalous dimension of the soft function and the function are known 

exactly to 3-loop order. The function §f(L) is known exactly only to order a s . Hoang and 
Kluth [15] argued that at order a 2 the function Sf(L) must be a polynomial of at most 2nd 
order in L, i.e. 'sf(L) = cf + c| L L 2 . In this paper, we show that this Hoang-Kluth ansatz 
does not hold; 'sf(L) is much more complicated. Certain moments of s/(L) contribute to the 
coefficients of 5(t) and 5(p) in the thrust and heavy-jet mass distributions. These moments 
were fit numerically in [15] and [12] using numerical calculations of the singular behaviour of 
these distributions in full QCD with the program event 2. In this paper, we produce these 
moments analytically and find that they are in excellent agreement with the most accurate 
available numerical fit [12]. 

Any L dependence at large L in 3/(L) turns into large logarithmic behavior of the hemi- 
sphere mass distribution (i.e. ln(M^/ 'Mr)). Since all of the p dependence is in s M (L), these 
large logs are not determined by RG invariance and correspond to so-called "non-global logs" . 
Dasgupta and Salam calculated the non-global logs for the related left-hemisphere mass dis- 
tribution in full QCD [26] and found no non-global logs (up to order L 2 ) for the CpUfTp color 
structure and an L 2 term with coefficient — for the CfCa term. We show below that the 
asymptotic behavior of ^f(L) in the full soft function is indeed of the form —^-L 2 for the 
CrCa color structure. We also find that both this color structure and the Cpn^Tp one have 
additional non-global single logs. These are especially interesting because the soft function is 
symmetric in L — > —L, which seems to forbid a linear term. The linear term appears through 
a complicated analytic function involving polylogarithms which actually asymptotes to \L\. 

This paper is organized as follows. In section 2 we review the factorization formula for 
the hemisphere mass distribution and its thrust and heavy-jet mass projections. Section 3 
computes the soft function in dimensional regularization. The calculation is complicated, so 
the results are summarized separately in 3.3. Section 4 discusses the result and presents the 
renormalized result for the integrated soft function, which can be compared directly to the 
predictions from SCET. Section 5 gives the previously missing terms in the singular parts 
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of the 2-loop thrust and heavy jet mass distributions, and compares to previous numerical 
estimates. Section 6 gives the full integrated hemisphere soft function which is compared to 
previous conjectures. The asymptotic form of this distribution, which exhibits non-global logs, 
is discussed in Section 7. Section 8 has some comments on predicting higher order terms with 
non-Abelian exponentiation. Conclusions and implications are discussed in Section 9. 



2 Event Shapes and Factorization in SCET 

The hemisphere soft function appears in the factorization theorem for the hemisphere mass 
distribution. The hemispheres are defined with respect to the thrust axis. Thrust itself is 
defined by 

T = max ( E jl P r nl ) , (2) 



where the sum is over all momentum 3- vectors in the event. The thrust axis is the unit 3- 
vector n that maximizes the expression in parentheses. We then define the light-like 4-vectors 

= (1, n) and n M = (1, — n). In the dijet limit T — > 1 and it is therefore more convenient to 
define r = 1 — T as the thrust variable so that r is small in the dijet limit. 

Once the thrust axis is known, we divide the event into two hemispheres defined by the 
plane perpendicular to the thrust axis. We define P£ and P£ to be the 4- vector sum of all of the 
radiation going into each hemisphere and Ml = yPf and Mr = v/Pj| to be the hemisphere 
invariant masses. When both Ml and Mr are small compared to the center-of-mass energy, 
Q, the hemisphere mass distribution factorizes into [8] 



(t dMf dM| 



H(Q 2 ,p) j dk L dk R J(M 2 L -Qk L ,p)J(M 2 R -Qk R ,p)S(kL,k Rj p). (3) 



Here, <Jq is the tree level total cross section. H(Q 2 ,p) is the hard function which accounts for 
the matching between QCD and SCET. J(p 2 ) is the inclusive jet function which accounts for 
the matching between an effective field theory with soft and collinear modes to a theory with 
only soft modes. Finally, the object of interest, S(kL, k R , /i) is the hemisphere soft function, 
which is derived by integrating out the remaining soft modes. 

In the threshold limit (small hemisphere masses), the thrust axis aligns with the jet axis 
and thrust can be written as the sum of the two hemisphere masses, 



Heavy jet mass p is defined to be the larger of the two hemisphere masses, normalized to the 
center of mass energy Q, 

p = ^max(MlM 2 R ). (5) 

When p is small, both hemisphere masses are small and the event appears as two pencil-like, 
back to back jets. 
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The factorization formula can be used to calculate thrust and heavy jet mass in the dijet 
limit as integrals over the doubly differential hemisphere mass distribution. Explicitly, 

and 

P (7) 
The thrust distribution can be written so that it depends not on the full hemisphere soft 
function but on the thrust-soft function, defined as 

S T (k, fjL) = J dk L dk R S(k L , k R , n)5(k - k L - k R ) . (8) 

Since the thrust soft function is dimensionless and its /i dependence is determined by renor- 
malization group invariance, the k dependence is also completely known. Thus at each order 
in a s only one number, the constant part, is unknown. In contrast, for the heavy jet mass 
distribution, the full ki and k R dependence of the soft function is needed for the factoriza- 
tion theorem. In particular, for resummation to N 3 LL order, only one number is needed for 
thrust (the constant in the 2- loop thrust soft function), which has been fit numerically, but 
for heavy-jet mass a function is needed [12]. In this paper we compute both the number and 
the function. 



3 Calculation of the Soft Function 

The soft function is defined as 

S(k L , k R , $ = 1-^2 S ( k R ~ n ■ P s R )$(k L - n ■ P S L ) (0| Y n Y n \X S ) (X s \ Y%y{ |0> , (9) 

c x s 

where P^ ,R is the total momentum of the final state \X S ) in the left and right hemisphere, 
respectively. The Wilson lines Y n and Y n are defined by 

Y^(x) = Pexp ^ig J ds n ■ A s (ns + x)^ Y^ n (x) = Pexp (ig J ds n ■ A s (ns + x)^ , (10) 

where P denotes path ordering and A s = A a s T a (A s = A a s T a ) are gauge fields in the fundamen- 
tal (anti-fundamental) representation. The soft function can be factorized into a perturbative 
(partonic) part and non-perturbative part which has support of order Aqcd [13]. 

The authors of [15] observed that the form of the soft function is constrained by the 
non-Abelian exponentiation theorem and RG invariance, which puts constraints on powers 
of logarithms of fi. The theorem also restricts the Cp color structure in the soft function to 
be completely determined by the one- loop result. Beyond this, however, the soft function is 
unconstrained. The one-loop calculation was done in [9, 27]. The main result of this paper is 
the calculation of the perturbative part of the hemisphere soft function to order a 2 s . Since the 
order a 2 s color structure Cp is given in [15], we will only calculate the CfCa and the CpUfTp 
terms. 
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3.1 CfCa color structure 

The order a 2 calculation involves pure virtual graphs, pure real emission graphs, and inter- 
ference between the two. The pure virtual contributions to the soft function give scaleless 
integrals which convert IR divergences to UV divergences, and are not explicitly written. The 
diagrams needed to compute the pure real emission contributions are shown in Figs. 1-4, 
whereas the interference graphs between the order a s real and virtual emission amplitudes are 
shown in Fig. 5. 

The integrals corresponding to the two diagrams in Fig. 1 (and the twin of diagram A ob- 
tained by interchanging k and q) are given, in d = 4 — 2e dimensions using the MS scheme, 
by 



2c 

An J J (27r) d J {2Ti) d k-k+q-q- 



^ 2 C F C A \ f d d q f d d k 1 , , 

2(-4^ 4 ) Cf-^ / 77^7 / — Hk L ,k R ) (11) 



and 



Ib = -4<? 4 



^ e iB\ 2i r &d q r d d k 1 



An J J (2n) d J (2vr) d (fc+ + q + ){k~ + q~) 
C F C A \ ( 1 1 \ 2 / 1 



c ^^)[y¥ + 7F¥) +c AiP¥ + ^)) F{k ^ (12) 

where k~ = n ■ k, k + = n ■ k and F{kii &r) contains the S(q 2 ) and S(k 2 ) factors which put the 
emitted gluons on shell and the phase-space restrictions in the definition of the hemisphere 
soft function. Explicitly, F{k,L, &r) is given by 

F(k L ,k R ) = ±(-2m) 2 5(k 2 )5(q 2 ) 

x 0(AT - k + )Q{q + - q-)5{k + - k R )5{q- - k L ) 

+ e(k+ - k~)Q(q- - q + )5(k- - k L )6(q + - k R ) 
+ 9(fc- - /c + )6(g" - q + )5(k + + q + - k R )5(k L ) 

+ Q(q + - q~)Q(k + - k-)5(k~ + q~ - k L )5(k R )\ . (13) 

In each diagram, the momentum has been routed so that the 4-vectors k and q correspond to 
the final state gluons. The gluonic contribution to the C f Ca color factor will be symmetric in 
k and q due to the fact that the radiated gluons are identical particles. In the first diagram, 
a factor of two has been added since the integrand is unchanged after k q, whereas in the 
second diagram, the SCET Feynman rules for two gluon emission from a single soft Wilson 
line automatically account for k -B- q. The factor of 1/2! needed for averaging over k -H- q 
is in F{k,L, k R ). Not shown in Fig. 1 is the graph that corresponds to the complex conjugate 
of diagram B. This diagram gives the same integral as diagram B. Since we are interested in 
the C f Ca contribution, the linear combination of interest is I a + 21 R . Diagram A and its 
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Figure 1: Diagrams A and B contribute to both C\ and CfCa color factors. 




Figure 2: Diagrams C and D contribute to the CfCa color factor. 




Figure 3: Diagrams E and F contribute to the CfCa color factor. 
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G H I 




Figure 4: Diagrams G, H and I. These classes of diagrams contribute to integrals Ig,Ih, and 
1 1 when the self-energy graphs involve gluons or ghosts and they contribute to integrals Ig,Ih, 
and 1 1 when the self-energy graphs involve a fermion/anti-fermion pair. 



J K L 




Figure 5: Diagrams J and K account for the interference between the one-loop virtual emission 
amplitude and the single gluon emission amplitude. Diagram L is the contribution from charge 
renormalizat ion . 
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identical twin are self-conjugate and only contribute once because they represent the squares 
of tree-level Feynman diagrams. 

There are four classes of diagrams involving the triple gauge coupling. Diagrams C and D, 
shown in Fig. 2, give the following integrals, 

I c = - g *C A C F (flfl) 2e I I d d k F(k L ,k R ) / Ar + 2<r + <T + 2t 



Air J J (2ir) d J (2ir) d (k~ + q~)(k + q) 2 \ k+q- q+k~ 

(14) 

/iV E \ 2e f d d q f d d k F(k L ,k R ) [2k + + q+ 2q + + fc+' 



Air J J {2ir) d J (2ir) d (k+ 
whereas, diagrams E and F, shown in Fig 3, give 

Ie = q^CfCa 



(15) 



/i 2 e 7B \ 2£ f d d q f d d k 



Air J J (2ir) d J (2ir) d 

1 1 \ q--k~ 



q- k~ J (k+ + q + )(k~ + q~)(k + q) 2 



F(k L , k R ) 



If — g A CpCA 



/x 2 e^ x 
Air 



d d k 



d d q 
(2ir) d J (2tt) 



1 

k+ 



k- 



(k+ + q+)(k~ + q-)(k + q) 2 



F(k L , k R ) 



Each of these diagrams has a complex conjugate and so they contribute twice. 

There are three self-energy topologies, shown in Fig. 4. The gluon and ghost self-energy 
graphs contribute to integrals Iq, Ih and h below. 



Ig = 9 A CfCa 



/i 2 e^ x 



d d q 



d d k 



1 



Air J J (2ir) d J (2ir) d (k~ + q~)(k+ + q+)(k + g) 4 
q + [(d - G)q~ -{d + 2)k~] + k + [(d - 6) AT -{d+ 2)q~] + 16k ■ q F(k L , k R ) 



Ih = q^CfCa 



u 2 e lE 
Air 



2c 



d d q f d d k 1 
j2nYJ (2ir) d (k- + q-) 2 (k + qY 



2{d + 2)q~k- -(d- 6)(k~) 2 -{d- 6)(g") 2 F(k L , k R 
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,> 2 e^V £ f d d g f d d k 

ii = g ^f^a 



An J J (27r) d J (2vr) d (k+ + q+) 2 (k + g) 4 

2(d + 2)g+£; + - (d - 6)(k + ) 2 - (d - 6)(g + ) 2 ] F(k L , k R ) 

As usual, cutting Feynman diagrams removes any symmetry factors that were associated to the 
cut lines prior to cutting. It is also worth reminding the reader that, in order to consistently 
combine the ghost emission diagrams with the gluon emission diagrams, we have to double- 
count the ghosts (they do not have the 1/2! symmetry factor that the gluons do). Diagram 
G has a complex conjugate graph which must be included but diagrams H and I, like diagram 
A, represent squares of tree-level Feynman diagrams and are therefore self-conjugate. 
Adding all of these contributions together, we have 

Sg A (k L , k R ) = I A + I H + h + 2(I B + I c + I D + I E + I F + I G ) (C F C A part) 

4r r /Ve^\ 2e / d d q f d d k f 2 

g^F^A^ ^ J J (2?r)(i J (2n) d \ (k ■ q) 2 k~ k+ q~ q+ (k~ + q~)(k+ + q+) 

x —k ■ q((k-) 2 q + (2k + + q + ) + 2k~q~ {(k + f - k + q + + (q + ) 2 ) + k + {q~f{k + + 2g + )) 
+ 2(k ■ qf(k~{2k + + q+) + q~(k + + 2g+)) 

+ (e - 1) ITif???fc?K»«'. (16) 

Before presenting the result for the C f Ca color factor, we briefly describe our general 
computational strategy. Normally, one expects scaleless integrals to be simpler than single 
scale integrals. In this particular case, the single scale integrals (with scale /cl/^j?) are actually 
much less technically demanding. This is true primarily because these contributions (see the 
first two terms of Eq. (13)) are integrable at e = 0. It turns out that this special feature of 
the problem more than makes up for the fact that single scale integrals are generically harder 
to evaluate than scaleless integrals. 

The calculation proceeds as follows for a single scale integral. First there will be an integral 
over angles (the integrand depends non-trivially on k ■ q) that can be done analytically to 
all orders in e. It is then convenient to Taylor series expand the resulting hypergeometric 
functions using the Hyp Exp package [28] for MATHEMATICA. In fact, the whole integrand 
can be expanded in a Taylor series in e and integrated term-by-term, due to the fact that 
the integral converges at e = 0. With a modest amount of knowledge of the basic functional 
identities satisfied by the polylogarithm functions, it is possible to do the resulting two-fold 
one-parameter integral in Mathematica and express the final result in terms of a minimal 
basis of transcendental functions. The results of our single scale calculations for both non- 
trivial color factors are tabulated in the Appendix. 

The evaluation of a scaleless integral (originating from the last two terms of Eq. (13)) 
begins in much the same way. Unfortunately, it quickly becomes clear that what remains after 
integrating over all angles has a non-trivial analytical structure (considered as a function of 
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e). In particular, the integral diverges at e = 0. Expanding an integral of this class under the 
integral sign is significantly more complicated and requires new tools. To begin, one should 
transform all hypergeometric functions in the integrand and expose their singularity structure. 
In this fashion, one learns that there is a line of singularities within the region of integration. A 
well-known procedure called sector decomposition [29] allows one to move singularities within 
the region of integration to singularities on the boundaries of the region of integration. Sector 
decomposition works as follows. Through a sequence of variable changes and interchanges of 
integration orders, all phase-space singularities are put into a canonical form. At this point, 
one can use an expansion in distributions to extract singularities in e under the integral sign. 
Finally, the entire integrand can be expanded in e in terms of distributions and ordinary 
functions and one can integrate the Laurent series term-by-term. 

Once we understood the computational procedure described above, it was straightforward 
to evaluate the integrals of interest for the CfCa color factor. The result has the form 

The first term corresponds to the first two terms in Eq. (13), those that account for the 
possibility that exactly one gluon is radiated into each hemisphere. It depends on fc A ( r , e), a 
dimensionless function of r = k^jk R and e. It can be written as an expansion in e as 

fo A (r, c) = f§l(r) + e/S(r) + e 2 /g(r) + • • • . (17) 

The expressions for f^ 1 '{r) are quite lengthy and are given in the appendix for n = 0,1,2. 
The second term in Eq. (17) accounts for the fact that both gluons can propagate into the 
same hemisphere and it has no non-trivial ki or k R dependence. gc A i e ) * s simply a constant 
with e expansion 

4 22 1 /134 4tt 2 \ II6C3 Utt 2 772 



'484C 3 4784 67vr 2 137vr 4 
9 + ~8T + ~27 90~~ 



The interference between the one-loop and tree-level single gluon emission amplitudes is 
shown in diagrams J and K of Fig. 5. The integrals associated with diagram J are scaleless 
and are set to zero in dimensional regularization. Diagram K gives the integral 



! K 



A(-g 4 )C A C 



d d q 1 f d d k 2q- - k~ 



(27r) d q-J (27r) d k+(q- - k~)(q - kfk 2 
x (-27ri)5(q 2 )[e(q- - q + )5(q + - k R )5(k L ) + 6(g + - q-)S{q~ - k L )5(k R )}. (19) 

There are 2 diagrams with the topology of diagram K. When they are considered with single 
real emission phase-space cuts, they can easily be mapped into each other and therefore give 
identical results. Both diagrams also have a complex conjugate graph and these obviously 
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give equal contributions as well. This is why Ix has an overall factor of 4 out front. After 
evaluating this integral, the real-virtual interference contribution becomes 



\ L Ft * 

where vc A (e) can be expanded in e as 



«c», (20) 



4 i 2tt 2 32C 3 tt 4 
: 30 



^(e) = ~ + ^ + ^-6^. (21) 



It is worth noting that, in this case, the application of the optical theorem for Feynman dia- 
grams is a bit subtle; one finds an explicit factor of exp(±27re) after doing the k integral (the 
sign of the phase depends on the precise pole prescription). Cutkosky's rules still apply pro- 
vided that one keeps only the appropriate projection of the complex phase. After a moment's 
thought it becomes clear that the real part, cos(7re) [independent of the pole prescription), is 
what one needs to keep to complete the calculation and derive the above result. 
The result of diagram L, including the complex conjugate graph, is given by 

S-feW = - {^)C F + J^-f, (22) 

where (3q = yC^ — f^/Tp is the first expansion coefficient of the QCD ^-function, (3(g) jg = 
f^Po- Finally, the total contribution to the C f Ca color factor is given by 

Sc A fa, k R ) = S« A (k Ll k R ) + S% A {h L , k R ) + S^(k L , k R ), (23) 
where S% e A n is the C F C A part of S Rcn . 

3.2 CftifTf color structure 

The diagrams involving a fermion loop contribute to the C pn fTp color factor and give integrals 
Ig, Ihi an d //• The first topology in Fig. 4, where the blob now represents a fermion loop, 
gives 



^ 2 e lE \ 2t f d d q f d d k 4(k + q- + k~q + - 2k ■ q) 



The phase-space cut is accounted for by 

F nf (k L ,k R ) = (-2m) 2 6(k 2 )6(q 2 ) 
x e(AT - k + )Q{q + - q~)5{k + - k R )5{q~ - k L ) + Q{k+ - k~)Q{q- - q + )6{k~ - k L )5{q + - k R ) 

+ e(AT -k + )Q{q~ -q + )5(k + + q + -k R )5(k L ) + Q(q + -q~)Q(k + -k~)5(k~ + q~ -k L )5(k R ) 
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The complex conjugate of this diagram gives the same result, so I g contributes twice. For the 
second and third topologies shown in Fig. 4, we get 

I H = gC F n f T F ^—) J— J J^ d {k+ + q+nk + y n f (k L ,k R ) (25) 

and 

'/i 2 e 7B \ 2£ f d d q f d d k -8k~q- 



Il = g * CF n f T F[ ^-) y^y___^_ F „ ;fe , M . (26) 

The sum of these contributions is 

SK(k L ,k R )=2I G +I H + l! 



g 4 C F n f T F 



fi 2 e^ E \ 2e f d d q f d d k 



An J J (27r) d J (2ii) d 



8_ / k + q + k q + - 2k ■ q _ k q _ k + q + 

{k + qY V (k+ + q+)(k- + q-) (k~ + q~y (k+ + g+) 2 ) n > 1 L ' R) 

Evaluating this integral gives 



S*(k L ,k R )= [^fc F n f T F 



As in the C f Ca case, the first term corresponds to the quark and anti-quark propagating into 
different hemispheres and it depends on r = ki/k R in a non-trivial way through a function 
fn f {r, e) • fn f (r, e) can be expanded in a Taylor series in e as 

/n / (r,e) = /(° / )(r) + e/( 1 / )(r) + ---. (28) 

The expressions for ft\ r) are given in the appendix for n = 0, 1. For the C F rifT F color factor 
n = 2 plays no role due to the fact that f!y(0) = fn/(oo) = 0. 

The second term Eq. (28) is present because both the quark and anti-quark may propagate 
into the same hemisphere as well. As before, this contribution has no non-trivial k^ or k R 
dependence. The constant g n , has a series expansion 

8 40 152 4tt 2 / 952 20tt 2 176Cs\ {oq . 

The final contribution to this color factor is from the charge renormalization, diagram L, 
the results of which were given in Eq. (22). Adding this contribution to the real emission 
contributions yields the final result for the C F rifT F color factor. It is 

S nf (k L , k R ) = S*(k L , k R ) + S^(k L , k R ), (30) 
where S^ n is the C F nfT F part of S Ren . 
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3.3 Summary of the Calculation 



In summary, we found that the 2-loop hemisphere soft function in d = 4 — 2e dimensions has 
the form 



S{k L ,k R ,fi) = (j^j 



k\ +2e v n> k R +2e K L> ) e 2 r(i 



(31) 



Here f(r,e) = /(1/r, e) is the opposite-direction contribution (where the two gluons or two 
quarks go into opposite hemispheres) and h(e) is the same-direction contribution. Since all 
the fi dependence is shown explicitly, h(e) cannot depend on k^ or k R by dimensional analysis. 
The second line is the contribution that comes from the interference of the first non-trivial 
term in the expansion of the charge renormalization constant and the 0(a s ) hemisphere soft 
function. It is proportional to fio = yCvi — ^T F Uf. 

There are 3 color structures, Cp,CpC 'a and CpnfTp. The Cp color structure is trivial - 
by non-Abelian exponentiation it is the square of the one-loop result. For the other two color 
structures the function f(r, e) is complicated. In both cases it is finite at e = 0, and in the 
CpUfTp case, / n .(0,e) = / n (oo,e) = 0. We write 

/(r,6) = /(°)(r) + e/( 1 )(r) + 6 2 / (2) W (32) 

The expansions in e of f(r, e) for the two color structures are given in the Appendix. Due to 
the fact that / re/ (0,e) = / n/ (oo,e) = 0, fn f (r) does not contribute to the renormalized soft 
function and is not given. 

For the same direction contribution, h(e), there are contributions from the real-emission 
diagrams and, for the C f Ca color structure, interference between tree-level real emission and 
one-loop real- virtual graphs. The real emission contributions we called g{e), and are given in 
Eqs. (18) and (29). The interference graphs are given by vc A (e) in Eq. (21). Adding these 
terms we get for the C f Ca color structure 

, 22 + OOA llvr 2 772 /484C 3 4784 67vr 2 14tt 4 



and for completeness, copying Eq. (29) 



40 152 4tt 2 / 952 20vr 2 176C 3 » 



4 Integrating the soft function 

Now we would like to expand and renormalize the soft function. At one-loop, all that is 
necessary for the expansion is the relation 



"1" 




In* 




-2e 




k_ 




k 







+ 



(35) 
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where the ^-distributions are defined, for example, in [9]. Unfortunately, this expansion cannot 
be used separately for k^ and k R , since the region where they both go to zero is not well-defined. 
For example, what does 8{kL)8{k R ) f {k L / k R ) mean? If we take ki — > first, then k R — > 0, 
then we pick up /(0). If we take /cl, k R — > holding k^ = k R , then we pick up /(l). Unless 
/(r) is constant, one must do the expansion more carefully. 

A simple solution is just to expand in distributions of p = kik R and r = k^jk R . This 
expansion is well-defined, and can be used to integrate any observable, such as thrust or 
heavy jet mass against the hemisphere soft function. For example, consider the integrated soft 
function: 

K(X,Y,fi)= [ dk L f dk R s(k L ,k R ,n). (36) 
Jo Jo 

This function contains the entire soft contribution to the integrated doubly differential hemi- 
sphere mass distribution. Since it is a function, rather than a distribution, we can use this 
integrated form to check the //-dependence and compare to previous predictions. 

We can calculate lZ(X,Y,ji) using Eq. (31) and the expansion in Eq. (35). For the same 
direction contribution (the real emission graphs, real/virtual interference graphs, and charge 
renormalization), the soft function is trivial to integrate in d dimensions. For the opposite 
direction contribution, the integral of the distributions is complicated by the overlapping singu- 
larities. It is a straightforward exercise in sector decomposition [29] to isolate the singularities 
and perform the integrations. The result can then be renormalized in MS. We find, for the 
opposite direction contribution, 



^(X,F, / ,) = (^) 2 |l / ( 2 )(0)-i/«(0)ln 



i 

2 





"1" 


f {1 \z)+ [dz 


In z 


/ dz 




to 


z 


+ Jo 


z 



f {0 \z)+\n 



XY 





"1" 


[ dz 
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z 



f (0 \z) 



Y/X 



dy 



Y/X dx /!%Kg (Q) 1 



xy 



(37) 



where /^(r) refer to the coefficients in the expansion in Eq. (32). The final compiled results 
for 1Z(X, Y, fj,) for the different color structures, including the same-direction and opposite 
direction contributions, are given in Sec. 6. 

The integrated soft function directly gives us the a 2 s soft function contribution to the 
integrated order o? s heavy jet mass distribution, 



da 



7 dp' = H{pQ,pQ,fi). 



(3* 



For thrust, the integrated distribution is not given in trivial way from the integrated soft 
function. However, it differs from the heavy-jet mass distribution only by a single finite 
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integral 



/<»>(./,) (39) 



xy 



which we can now compute for the C f Ca and C F njT F color structures. Adding also the C F 
terms, which were already konwn, the result is 



R t (t,h) = R p (t,h) + 



K -C 2 F -} ( ~ - 8C3 ) C F n f T F 



(40) 



45 r \3 

+ ( 32Li 4 l + 22C 3 + 28C 3 ln2 - | - ^ + ^ - ^ 2 ln 2 2 ] C F C, 

5 Numerical check for thrust and heavy jet mass 

As a check on our results, we can use the soft function to calculate the soft contribution to the 
differential thrust and heavy jet mass distributions. The singular parts of these distributions 
at 0(a 2 ) were previously determined up to four numbers: the coefficients of 5{t) and 5(p) 
for the CfCa and CprifTp color structures. Until now these four numbers were unknown 
and had to be fit numerically using the event 2 program. We can now use our results for 
the hemisphere soft function to replace these numerically fit numbers with analytical results. 
The coefficients of the ^-functions are the same as the constant terms in R p {p) and R t (t), for 
which formulae were given in the previous section. 

The unknown soft contributions to the coefficients of S(r) and 5(p) were denoted cf and 
4 p in [12]. We find 



/ 2032 871tt 2 16vr 4 4 In 4 2 4 , 

cf„ = —Cl + + + -7T 2 n 2 2 - 28C 3 n2 

2p 2 F l 81 54 9 3 3 



+ + (42) 

These numbers were fit numerically in [11, 12, 15] based on a method introduced in [11]. 
The procedure involves subtracting the singular parts of the thrust and heavy jet mass dis- 
tributions, which are known analytically from SCET, up to delta-function terms, from the 
full QCD distributions for thrust and heavy jet mass calculated numerically with the program 
event 2. The difference is then integrated over and compared to the total cross section, 
which is known analytically, minus the analytic integral over the singular terms. The highest 
precision fits were done in [12] so we compare only to those. The result is 

4 = (48.7045)C£ + (-56.4990)CVC A + (43.3905) C F n f T F (analytic result) 
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= (49.1)C£ + (-57.8)C F C A + (43A)C F n f T F (fit result [12]) (43) 

and 4 p 

cf p = (48.7045)C| + (-33.2286)C F C A + (50.3403) dm f T F (analytic result) 

= (49.1)C| + (-33.2)C F C A + (50.2)C F n f T F (fit result [12]). (44) 

The percent errors for these numbers are 0.8%, 2%, 0.02% for cf and 0.8%, 0.08% and 0.2% 
for cf p respectively, with an average error of around 0.5%. This is excellent agreement. Note 
that the Cp terms were already known when the fits were done, so small errors were expected. 

For completeness, the complete contributions of the soft function to S(p) and S(r), denoted 
by Dg and Dg at order a 2 s are 

^ - {* - \^ + ™ i^^-^y °«*> {-^ - ^) }■ 

of - {<, - + M (»& + 2£ - £) + «w (-^ - "£) }. 

(45) 

One can also use cf and cf to get the complete coefficient of 5(r) and S(p) including jet and 
hard function contributions, using Appendices C of Refs. [11] and [12]. 



6 Hemisphere mass distribution 

The numerical check performed in the previous section provides strong evidence that our 
analytical results are correct. With these results in hand, we can now compare to other 
features of the hemisphere mass distribution and the integrated hemisphere soft function, 
K(X,Y,p). 

The entire /x-dependence of 1Z(X, Y, p) is predicted by SCET. Indeed, renormalization 
group invariance predicts that the differential soft function must factorize in Laplace space, 
as in Eq. (1). The Laplace transform is defined by 

S(x L ,x R ,fi)= / dk L / dk R S{k L ,k R ,p)e-^ k ^ B e- XRk ^\ (46) 
Jo Jo 

where the je factors are added in the definition to avoid their appearance elsewhere. The 
factorization theorem then implies 

s(x L ,x R ,p) = s IM (\nx L p)s^(\nx R p)s f (x L ,x R ). (47) 

The RG-kernel s^(L) is determined by the renormalization group invariance of the factorization 
formula, and is expressible in terms of the anomalous dimensions of the hard and jet functions, 
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which are known up to a^. The finite part Sf(L), until now, has been known only to a s . This 
Laplace form leads to a simple expression for the integrated soft function in SCET [12] 



K(X,Y,fj) = s{d m ,d m ,n) (— 

V A 4 



-j Em 



>I2 



-1EV2 



nv2 



??l=»72=0 



(4£ 



The //-dependent terms in the order a 2 s integrated soft function calculated in this way agree 
exactly with the //-dependent terms in 1Z(X,Y, ji). In fact, it is helpful to separate out those 
terms. To that end, we write the a? terms as 



M - - 



k>4 



(49) 



where 1Z^(X/ \x,Y j \x) is the part coming directly from the s M (L) terms and TZf(X/Y) is the 
remainder, which comes from 'sf(xL,xn). The result for TZ^X/fi, Y/fi) is 



X Y 
H ' fJL 



X 20 
fx 3 
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/' 



; l n 4 11 _ ^ ln 2 H + 16 ln 2 H ln 2 



(>1 xy o , 4 y 20 2 , 2 y 28vr 4 
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ln d 



+ -7r 2 ln 2 



X 26£ 



-28C 3 ln 



/t 3 
XY 8 
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X 22 
7T ~ ~9~ 



! 2l XY 808, XY 
-it In — — H — In — — 



27 



, , Y 4 2l 2 y 268, 2 y 352C 3 4tt 4 268tt 2 
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CfCa 



+ 



32 , oX 80, 2 X 8 9 , XY 224, Xy 
ln d h — In 2 h -vr 2 In — — In ■ 



fi 9 /i 9 



27 



807T 2 



32, oY 80, 2 Y 128C 3 

— In 3 — + — In 2 — 

9 ijl 9 fi 9 27 



CpTpUf. 



(50) 



The part of the soft function not determined by RG-invariance is represented entirely 
by t>f(xL,Xii). This function is //-independent and can only depend on the ratio Xl/xr by 
dimensional analysis. Moreover, it is symmetric inx^ Xr, since the hemisphere soft function 
is symmetric in k^ <H- k^. Hoang and Kluth claimed [15] that it should only have logarithms, 
and up to order a 2 , only have In and In 2 terms. Their ansatz was that 



Sf(x L ,X R ) 



Hoang-Kluth 



Ait 



(51) 



with cf = —Cftt 2 already known. 
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To check the Hoang-Kluth ansatz, the easiest approach is to look at the contribution of 
Sf{xL,XR) to TZ(X,Y, fi), which we called lZf(X/Y). For the Hoang-Kluth ansatz, the result 
is 



^(^Hoang-Kluth = Q S + ^(^2 z _1L), 

3 



(52) 



The values of cf and cf L which get right the singular parts of the thrust and heavy jet mass 



distributions are given in Eqs. (41) and (42) with cf L = J?(cf — cf 



The exact answer, at order a 2 s is 
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3 3 3 9(2 + 1) 



506C(3) 16tt 4 871tt 2 2032 
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54 



32 

32Li 3 (-z)-— Li 2 (-z)ln(z) 



8 2-I In (z) 16 2l , N 184(3 154tt 2 136 

+^—, — r H 7T In U H H 

3(^ + 1) 9 V ; 9 27 81 

This is clearly very different from the Hoang-Kluth form. 



(53) 



7 Asymptotic behavior and non-global logs 

The factorization theorem is valid in the dijet limit when the hemisphere masses are small 
compared to Q; however, there is no restriction on the relative size of the two masses. In 
addition to logarithms In Ml,r required by RG invariance, there may be logarithms of the 

form In that enter at order a 2 . These logarithms cannot be predicted by RG invariance 
and are known as non-global logarithms. Salam and Dasgupta have shown that non-global logs 
appear in distributions such as the light jet mass. They argued that in the strongly-ordered 
soft limit, when M L <C M R < Q, the leading non-global log should be ~(j^) 2 ^f- C F C A In 2 
in full QCD. This double log was reproduced in [30]. 

Non- global logs must be present in SCET, since for small Ml and Mr, the entire distri- 
bution is determined by soft and collinear degrees of freedom. The non-global logs cannot 
come from the hard function, which has no knowledge of either mass, or the jet function, since 
each jet function knows about only one mass. Thus, they must come from the soft function. 
Moreover since, by definition, they are not determined by RG invariance, they must be present 
in the /i-independent part, 1Zf(X/Y) of the integrated hemisphere soft function, TZ(X,Y, fj,). 
This function was given explicitly in Eq. (53). 

To see the non-global logs in IZf(z) we can simply take the limit z — > 00. Note that 
TZf{z) = TZf(-) so this is also the limit z — > 0. The asymptotic limit of TZf(z) for large or 
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Figure 6: The contribution of the part of the soft function not fixed by renormalization group 
invariance to the hemisphere mass distribution, TZf(z) is shown. On the left is the CpUfTp 
color factor and on the right is CfCa, both as a function of In z = In y. The solid black curve 
is the exact result of Eq. (53). The dashed red curve is the plot of the small lnz expression of 
Eq. (55) and the dotted blue curve give the large In z behavior of Eq. (54). The kink is due 
to a sign flip, since the linear term appears as | hiz\. 
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(54) 
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CfCa- 



There are two important features to note in this expansion. First of all, in the CfCa color 
structure there is a term — In 2 z, which is the leading non-global log found by Dasgupta 
and Salam and in [30]. But we also see that there are sub-leading non-global logs, of the 
form | In z \ . The absolute value is necessary to keep the expression symmetric in z — > - . It is 
interesting to see how this sign flip comes out of the full analytic expression. 
Next, let us look at z ~ 1. Here we find 
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(55) 



We see there is a double logarithmic term for both the CfCa and CpfifTp color structures. 
This is consistent with an analysis performed in [15] of an observable p a = max(aM|, M^)/Q 2 . 
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They found that the integrated p a distribution looked like In 2 a for a ~ 1. This quadratic 
behavior in In a corresponds exactly to the quadratic behavior in the z ~ 1 limit in (55). 

We show in Figure 6 the exact finite function TZf{z) and a comparison to its asymptotic 
behavior at small \nz and large lnz for the CfCa and CpnjTp color structures. For both 
color structures, the exact curve is well approximated by a parabola for small In z. At large 
lnz, for the CpUfTp color factor, the exact result approaches a linear function whereas the 
CfCa color structure has In 2 z dependence with a different coefficient than for the small lnz 
limit. The CfCa term has a linear term as well. 

As we have discussed, the integrated hemisphere soft function contributes directly to the 
doubly differential hemisphere mass distribution. In the limit where both hemisphere masses 
are small, and well separated, the soft function gives the dominant contribution. In this 
regime, we can read off that the leading non-global logarithms are given by 7?.j >1 (M|/M|,) 
in Eq.(54). The In 2 term has an identical coefficient to that found in [26]. The subleading 
non-global logarithm is a new result. 

8 Exponentiation 

These non-global logarithms become important when one scale becomes parametrically larger 
than the other. The separation of scales suggest that at the higher of the two scales, one 
may be able to match onto a new effective theory and then run the matching coefficient 
between the two scales. In fact, the — In 2 z term in this calculation has its origin in /(0, 0), 
where f{r,e) is the opposite-direction contribution to the hemisphere soft function, as in the 
Appendix. Indeed, we found that 

fc A (0, e) = ^ + O(e), / B/ (0, e) = G(e 2 ), (56) 

which is consistent with the leading non-global logarithm only having the CfCa color struc- 
ture. Since it is the e° part of this expression which contributes, and there are double soft 
poles, the full expansion also has terms like /(0, 0) In 2 Thus /(0,0) can be thought of as 
an anomalous dimension, providing hope that these non-global logs might be resummed in an 
effective theory. A consistent framework may require some kind of refactorization, like the one 
found for a related event shape, r w , in [25]. Ideas along these lines were suggested in talks by 
Chris Lee [30, 31]. Lee and collaborators proposed that the leading non-global logs might be 
resummed with effective field theory, although no details were given. 

There is actually good reason to believe the resummation of non-global logs is more chal- 
lenging than the types of resummation done in SCET. To see this, we first consider the 
predictions from non-Abelian exponentiation. Non-Abelian exponentiation applies only to 
the case of pure QCD, without fermion loops. In this case, it says that the full soft function, 
in Laplace space, can be written as an exponential of 2-particle irreducible diagrams. At order 
a™, new contributions can appear only to the maximally non-Abelian color structure, CpC 1 )^ 1 . 
For example, at two-loops, this tells us that the Cp color structure is given entirely by the 
exponential of the one-loop Cf color structure. At 3-loops it predicts the entire Cp and C f Ca 
color structures. 
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To be more specific, the soft function in Laplace space factorizes as in Eq. (47), with the 
terms and the s/ terms separately exponentiating, as explained in [15]. So we can write 



8f[x L ,x R ) = exp 



a 



7 J' nf (xl,xr) 



+ C F C A s 



(57) 

where lfj' n} ' and sf' Cj ^ are the Laplace transforms of the Cpn^Tp and CfCa color structures 
in the 2- loop soft function. 1 Such a rewriting has no content unless there is some restriction 
on the terms appearing in the exponent. Non-Abelian exponentiation tells us that the higher- 
order terms with Cf and C A s only must be maximally non-Abelian, CfC^~ . 

This implies, for example, that at 3-loops we know 2 color structures. Explicitly, 
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There are 4 remaining color structures, CfC a , Cpn 2 T F , CpC A^fTp and CpnjTp which are 
still unknown. Actually, the Cpn 2 T F color structure at 3-loops should not be hard to compute, 
but there is no known general formula for how the rif color structures exponentiate (see [35] 
for some discussion). 

From the exponentiation formula, one can read off the missing parts of the soft contribution 
to the 3-loop thrust and heavy-jet mass distributions. Indeed, for n > 2, we have 
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(59) 



and similarly for c 



2 

npi 



with cf and 



2/< 



given in Eqs. (41) and (42). These constants can be 



included in future a s fits or, once the finite part of the 3-loop jet function is computed, 
compared to extractions from the full thrust distrubtion at NNLO [36]. 

Returning to the exponentiation of non-global logs, recall that the leading non-global log 
comes from s^' Ca \xl,xr) = —^-L 2 , with L = hi§|. Thus non-Abelian exponentiation 
predicts a series with terms (C 'pC AOi 2 L 2 ) n ~ l , as well as cross-terms with the Cf one-loop color 
structure which are subleading. The question is whether this is the entire resummation of 
the leading non-global log. It seems like the answer is no, since there is no apparent reason 
why 3-loop graphs cannot produce terms which scale like o? s CfC 2 a Ij Z (or even a^L 4 ) for large 
L. A clue that these terms do exist comes from the numerical resummation of the leading 
non-global log at large -/V c in [26]. These authors found that the resummed distribution could 
be fit by an exponential, but it is numerically different from the pure {C pC AOi 2 L 2 ) n ~ 1 terms 
predicted by Eq. (57). Since Cp and Ca both scale as N c at large N c , this implies that there 
must be a o? s CfC 2 a I? term at 3-loops (and no a^L 4 term). Thus the resummation of even the 
leading non-global log may require a way to predict arbitrarily complicated color structures. 
It would be exciting to see how this can be done in the effective field theory framework. 



1 Although we have not computed the Laplace-space soft function directly, it was calculated by another 
group after the first version of this paper appeared [34]. It has a qualitatively similar form to TZf(^). 



21 



9 Conclusions 



In this paper, we have presented the complete calculation of the hemisphere soft function to 
order a 2 . This is the first 2-loop calculation of a soft function which depends on two scales 
in addition to the renormalization group scale \i. The hemisphere soft function, S(Icl, /i), 
depends on the components of the momenta going into the left and right hemispheres. In a 
one-scale soft function, such as the Drell-Yan soft function, S , £ ) y(fc,/x) [22, 23, 32], the thrust 
soft function SV(£;,/i) [9, 11] or the direct photon soft function S^(k,fi) [16], all of the k 
dependence is fixed once the /i-dependence is known. Since the //-dependence is fixed by RG 
invariance, these functions are often completely determined. For multi-scale soft functions, 
like the hemisphere soft function, there can be additional dependence on the ratio r = k^fk^. 
We worked out this dependence explicitly at order a 2 , and the result is more complicated than 
previously anticipated. 

We performed a number of checks on our calculation. The /i-dependence of the result was 
entirely known by virtue of the factorization theorem in SCET, and we have confirmed that the 
/i-dependence of our hemisphere soft function matches the result obtained from factorization 
analysis. In addition, the result allows us to produce analytic expressions for all of the singular 
terms in the 2-loop thrust and heavy jet mass distributions. The constant terms in the singular 
distributions were previously unknown and had to be extracted from numerical fits [11, 15, 12]. 
We found our analytical results to be in excellent agreement with the very precise recent 
numerical fit of [12]. 

The full hemisphere soft function produces the leading and sub-leading non-global logs in 
the hemisphere mass distribution. Previously, only the leading double-log term was known, 
from a calculation in the soft limit of full QCD [26]. In this work we reproduced that double 
logarithm and, furthermore, showed the existence of a sub-leading single logarithm. This single 
logarithm, of Ml/ Mr, is interesting because In Ml/ Mr seems like it should be forbidden by 
the Ml -B- Mr symmetry. Curiously, we find that the complicated behavior of the hemisphere 
mass distribution when Ml ~ Mr allows the single log to flip sign and it manifests itself 
as ln[max(M^, Mr)/ min(ML, Mr)} = \IuMl/Mr\. Our calculation is the first to exhibit a 
sub-leading non-global logarithm of this type. 

Besides being of formal interest, the hemisphere soft function at 0(a 2 ) is a crucial compo- 
nent of the resummed heavy-jet mass distribution at N 3 LL order. Previous fits to a s at this 
order assumed a simple form for the soft function, using the Hoang-Kluth ansatz. We have 
shown that this ansatz is valid only in the limit that kL ~ kR. With the exact 0(a 2 ) soft 
function in hand, one source of uncertainty in the a s fits to event shapes can be removed. 

This work also has implications for calculations of distributions at hadron colliders. At 
hadron colliders, there are necessarily many more scales in relevant observables than at e + e~ 
machines. For example, jet sizes and veto scales play a critical role in many analysis [33, 19, 
25]. For multi-scale observables to be computed in effective field theory, we need a better 
understanding of multi-scale soft functions, such as this exact result on the 2-loop hemisphere 
soft function provides. 
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A Opposite direction contributions 

The following are the first 3 terms in the e expansion of fc A from Eq. (17). 

r(0) , _ /r(llr 2 + 21r + 12)ln(r) 7r 2 (r + l) 2 + 2r 2 



fr(r) = 8 — -, ' ' — + v ' ' + ln 2 (r + 1) 

JCaK ! V 3(r + l) 3 3(r + l) 2 V ; 

-ln(r)ln(r + l)-yln(r + l)j , (60) 



„m, , 8(-llr 3 -9r 2 + 9r + ll)Li 2 (-r) , . , , _ 

f [ cl(r) = ~ 3( r + i)3 + 24Li 3 (-r) - 16Li 2 (-r) log(r) 

4r (llr 2 + 21r + 12) log 2 (r) 8r (67r 2 + Ulr + 60) log(r) 



3(r + l) 3 9(r + l) 3 

4 (r 3 (IItt 2 - 36C(3)) + r 2 (-108C(3) + 32 + 21tt 2 ) + 4r (-27C(3) + 8 + 3tt 2 ) - 36((3)) 

9(r + l) 3 

8(-llr 3 -9r 2 + 9r + ll)ln(r + l)ln(r) A _. lNl2/ , 4 . o2 m/ 

+ — -, tt^ — — — - 4ln(r + 1) ln 2 (r) + - (134 + 3vr 2 ) ln(r + 1), 

o(r + iy 9 

(61) 

(2) _ 8 (67r 3 + 81r 2 - 81r - 67) Li 2 (-r) _ 8 (55r 3 + 117r 2 + 81r + 11) Li 3 (-r) 
JCa^) - 9(r + l) 3 3(r+ l) 3 

32 (llr 3 + 9r 2 - 9r - 11) Li 3 (^) 16 (llr 3 + 9r 2 - 9r - 11) Li 2 (-r) ln(r + 1) 
3(r + l) 3 3(r + l) 3 

8(33r 3 + 75r 2 + 57r + ll)Li 2 (-r)ln(r) m f 1 *\ f r 

16Li 4 I 7- I — I6L1 



3(r + l) 3 \r + lj \r + l 
- 8Li 2 (-r) ln 2 (r) + 16Li 2 (-r) ln(r) ln(r + 1) + 8Li 3 (-r) ln(r) + 16Li 3 f ) ln(r 



+ 



x r + 1 

4 (-40 (4r 2 (27C(3) - 2) + r(189C(3) - 8) + 99C(3)) + 5yr 2 r (67r 2 + 147r + 66) + 33vr 4 (r + l) 3 ) 

135(r + l) 3 
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32(12r 2 + 21r + ll)ln 3 (r + 1) 4r (llr 2 + 21r + 12) ln 3 (r) 4 



+ 



9(r+l) 3 9(r + l) 3 3 

16 (12r 2 + 21r + 11) ln(r) ln 2 (r + 1) 4r (67r 2 + Ulr + 60) ln 2 (r 



hr(r) ln(r + 1) 



+ 



3(r+l) 3 9(r+l) 3 

16r (193r 2 + 384r + 177) ln(r) 4 (-llr 3 - 9r 2 + 9r + 11) ln 2 (r) ln(r + 1) 
27(r + l) 3 + 3(r + l) 3 

8 (tt 2 (66r 3 + 90r 2 + 9r - 33) + 2 (193r 3 + 561r 2 + 561r + 193)) ln(r + 1) 

27(r + l) 3 

8(67r 3 + 69r 2 -93r + 37r 2 (r + l) 3 -67)ln(r)ln(r + l) 4 l4 . . 

H ; In \T + 1) 

9(r + l) 3 3 V ; 

32r ln 2 (r + 1) 

+ 8 ln 2 (r) ln 2 (r + 1) + 16C(3) ln(r + 1) - 16C(3) log(r) + V ; . (62) 

3(r + 1) 

The following are the first 2 terms in the e expansion of f nf from Eq. (28). 

(0) 16 (2r(r + 1) - 2(r + l) 3 ln(r + 1) + r(r(2r + 3) + 3) ln(r)) 

" 3(r + l) 3 ' [b6) 

/W(r) = (-12 (r 3 - 1) Li 2 (-^ - 32r 3 ln(r + 1) + 3ttV + 20r 2 

-96r 2 ln(r + 1) - 3 (4r 3 + 3r 2 + 3r - 2) ln 2 (r) + 4 ln(r) (3 (r 3 - l) ln(r + 1) 
+r (8r 2 + 21r + 3) ) + 3?r 2 r + 20r - 96r ln(r + 1) - 32 ln(r + 1) + 2tt 2 ) . (64) 
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